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1^"^ 1 We provide a microscopic derivation for the non-Markovian master equation for an atom- 

' cavity system with cavity losses and show that they can induce population trapping 

QQ ' in the atomic excited state, when the environment outside the cavity has a non-flat 

f**^ , spectrum. Our results apply to hybrid solid state systems and can turn out to be helpful 

to find the most appropriate description of leakage in the recent developments of cavity 
quantum electrodynamics. 
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1. Introduction 

Cavity quantum electrodynamics (CQED) is one of the most important fields of 
quantum optics, both from a fundamental and an applicative point of view "'^. In- 
deed, the interaction of an atom and a single mode of the quantized electromagnetic 
field inside a high-Q cavity can be exploited to study the properties of highly non- 
classical states of the atom-cavitysystem and to use their non-classical features for 
quantum information processing^. In this context it is fundamental to take into 
account the unavoidable coupling between the atom-cavity system and the envi- 
ronment external to it. In general, the coupling with an environment gives rise to 
dissipation and decoherence phenomena , which can stro ngly affect the dynamics of 
the quantum system under study. In previous paper J^El ^e derived a master equa- 
tion for an atom-cavity system which takes into account from the very beginning a 
coupling between the atom and the cavity described by the Jaynes-Cummings (JC) 
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model^ . The approach followed in Refs. |4| and [5| is different from the one usually 
reported in the literature to describe cavity losses. Indeed the usual model of cav- 
ity losses has been introduced in a phenomcnological way, i.e., one microscopically 
derives the master equation for the cavity only, described as a quantum harmonic 
oscillator coupled to a bosonic bath, and then assumes that the presence of the 
atom inside the cavity does not affect the structure of the dissipator in the master 
equatioiPl. Through the microscopic model derived in Refs. IH |5] it is possible to 
give a complete justification of the phenomcnological model, which turns out to be 
valid when the spectrum of the environment is flat and its temperature is zero. The 
situation is quite different if the environment has a non-flat spectral density. This is 
not an uncommon situation nowadays: indeed there are new-generation CQED ex- 
periments which are performed on high-Q cavities created inside photonic bandgap 
(PBG) materialJ^. For this system we proposed a non-Markovian model of cavity 
losses^l which predicts a dynamics completely different from the one predicted by 
the phenomcnological model, giving rise to phenomena such as population trapping 
due to cavity losses. Scope of this paper is to provide some details of the derivation 
of the non-Markovian master equation for the JC model and to show that the main 
predictions of the theory, namely population trapping due to cavity losses^, does 
not depend on the shape of the spectral density one chooses, but only on the values 
of the asymptotic decay rates corresponding to the transitions of interest. We clar- 
ify this point by taking into account two different models of environmental spectral 
densities and showing that, for the same asymptotic values of the decay rates, the 
predictions differ only in the details of the very short-time dynamics, while the 
two densities give the same predictions for longer times. The paper is structured as 
follows. In Sec. 2 we present the non-Markovian model of cavity losses, in Sec. 3 we 
present the calculation of the time-dependent decay rate for a Lorentzian spectrum 
and summarize the dynamics of the system in this case. In Sec. 4 the case of a 
spectrum with a Lorentzian gap is presented along with some conclusive remarks. 

2. The non-Markovian master equation for cavity losses 

The system we study consists of a two-level atom interacting with a mode of a 
cavity coupled to a bosonic environment. Calling loq the Bohr frequency of the 
atom and \g) and |e) its ground and its excited states respectively, the interaction 
between the atom and the cavity mode is described, at resonance and in units of 
/i, by the JC HamiltoniarP Hjc = ^'^z -I- + ^ {'^'^+ -I- aV_) , where (a) 

is the creation (annihilation) operator of the mode, a_ — \g) (e|, a+ = \e) {g\, and 
a^^\e){e\-\g){g\. 

The cavity mode interacts with a bosonic reservoir through an interaction 
Hamiltonian linear in the bosonic operators of the reservoir and of the cavity mode. 
More precisely we assume that the dynamics of the total system (atom-cavity sys- 
tem and environment) is described by a Hamiltonian H = Hs + He + i^int, where 
Hs — Hjc is relative to the atom-cavity system. He = X]fe '^fe^I^fc is relative 
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to the environment and Hint — {o- + o,^) 9k y>k + b\j describes the system- 
environment interaction, with the frequencies of the environment osciUators, bj, 
(&fe) the creation (annihilation) operator of quanta in the k-th environmental mode, 
and gk the coupling constants. Since the reservoir causing cavity losses is immersed 
in the PBG material, we expect its spectrum to be non-flat: then, to be rigorous, 
the master equation must be derived in the framework of a non-Markovian theory. 

Using the second order of the time-convolutionless (TCL) expansiorPof an exact 
non-Markovian master equation, and neglecting the atomic spontaneous emission 
and the Lamb shifts, it is possible to show that for a weak enough damping the 
master equation for the system is local in time and that its structure, in the rotating 
wave approximation (RWA), is the same of the Markovian master equation for the 
same system, with the important difference that now the decay rates are time- 
dependent. Examples of non-Markovian master equations of this type can be found 
for example in Refs. [31 [TUlandfTTl In our analysis we will focus on the case of one 
initial excitation only and we will consider a reservoir at zero temperature. Then, 
the non-Markovian master equation for the atom-cavity system density operator p, 
in this case, is the following one—: 

pit) = -i[Hjc,p] 



+7 (wo + t) Q \Eo) | p{t) {Eo\ - i {E,^+\ , p{t)} 

+7(c^o - ^,t) Q \Eo) {E,.-\pit) \E,^.) {Eo\ - {E,,-\ , p{t)}^ , (1) 

where |-Ei,±) = ± |0, e))/-\/2 are the eigenstates of Hjc with one total ex- 

citation, with energy a;o/2 ± fi, and |-Eo) = |0,.g) is the ground state, with energy 
—ujo/2. The time-dependent decay rates for l-Ei,-) and are 'y{u!o — il,t) 

and 7(a;o + ^, t) respectively. Equation ([T]) is the non-Markovian extension of the 
Markovian model introduced, for the same case, in Ref. i4j 

By direct substitution, it is possible to show that, if the system starts from the 
state |0, e), i.e., if the atom is initially excited and the cavity is initially empty, the 
solution of Eq. ^ is the following one: 

/ 1 i-(t) 1 ^+(*)\ 1 i-W 

Pit) = 1 - -e-^ - -e-^ \Eo) {Eo\ + -e"^ \E^,.) {E,,.\ 



2 2 

1 i+(t} 1 -f^(t)+i"+(t) 



-e-^ \Er,+) {E,M - -e 1 (e^^" \E^,^) {E^M + h.c.) , (2) 

where I±{t) — /p 7(079 ±f^, t')dt' . From Eq. ^ one can compute all the populations 
we are going to show in the following. Below we will study the behavior of the non- 
Markovian time-dependent rates 7(079 ± which, through the quantities I±{t), 
lead to non-Markovian behavior and to population trapping. 
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3. The non-Markovian decay rates for a Lorentzian spectrum 

As a first model of environment at zero temperature witli non-flat spectrum, we 
consider the Lorentzian distribution^ 

1 aX^ 

J{^) = TTi \2~rT2' 

where a is the system-environment coupling strength, and A is the width of the 
distribution, describing also the inverse of the reservoir memory time. The case of 
Lorentzian spectrum is analytically treatable, while capturing important features 
of the non-Markovian dynamics we are interested in, i.e., the time-dependence of 
the decay rates and their different stationary values. We consider the case in which 
the spectrum is peaked on the frequency of the state i.e., loi = cjo — ^i 

where ujq is the atomic Bohr frequency and is the Rabi splitting due to the JC 
interaction. 

The rate 7(w,i) for a generic transition with Bohr frequency lo is equal to 
7(cj, t) = 2Re {r(u;, t)}, where T{lo, t) is related to the spectral density J(w) through 
the relatiorM 

pt r+oo 

T{LJ,t)= dr duj'S'^-'^'^^J{uj'). (4) 



Jo J-oo 

This relation gives the right Markovian decay rates and Lamb shifts for t — > +oo. 
By performing first the integral with respect to r, one obtains: 

^^^r{i^-i^')r ^ Sinjuj' - Uj)t _ , 1 - C0S(CJ^ - Lu)t 



(5) 



which, substituted into Eq. (|3]), gives: 



r{uj,t)^ duj'J{uj') ^ --i duj'J{uj') '-. (6) 

J-oo LU'-UJ UJ'-LU 

The second term in Eq. ([6]) gives a time-dependent Lamb shift which in the following 
will be neglected. The first term, i.e., the real part, is half the time-dependent decay 
rate 7(0;, t). 

Specializing to the Lorentzian spectral density given in Eq. ([3]), the real part of 
Eq. ([6]) becomes: 

Re {r(a;, t)} = -Im e^ * / dcu'- , ^ , (7) 

whose last form is suitable for an evaluation by means of the method of the residues, 
by closing the integration path with a circle of ray R on the upper complex half- 
plane, where the integral vanishes when R — > oa^. Such an evaluation is straight- 
forward. It is only worth noting that the pole on the real axis, at lu' — lu, must be 
circumvented from below, so that its residue is taken for half its value, with positive 
sign. The reason of this choice is that, as we will see, it gives a positive stationary 
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value for the decay rate, while the other choice would give a negative stationary 
decay rate, which would be unphysical. 

By taking twice the quantity in Eq. ([7]) and evaluating the integral, we finally 
obtain the following expression for the decay rate j{uj,t): 



■ sin(wi — uj)t — cos(cJi — uj)t 



.(8) 



From Eq. ([5]) we clearly see the general behavior of the time-dependent rates: all 
the rates ^{uj, t) are zero at t = 0, then they oscillate in time with a time-dependent 
average value, till they reach stationary values for t ^ A""'^. These stationary val- 
ues are proportional to J{ijj), i.e., they are equal to the rates one obtains from a 
Markovian theory. For these reasons, as anticipated, the quantity A~^ can be seen 
as the memory time of the system-reservoir interaction and non-Markovian effects 
are expected to occur for times shorter than A~^. 

In particular, taking the peak of the spectrum in uji = uja ^ ^ a-nd substituting 
w = Wo ± f^, one obtains the decay rates for the two dressed states |-Ei,±), which, 
along with the solution of the master equation in Eq. ^ , allow us to compute the 
dynamics of the atom cavity system presented in Ref. [9l to which we refer for the 
plots of the populations considered. The population of the ground state of the atom- 
cavity system \Q,g) increases quadratically for short times, with some oscillations 
superimposed which are signatures of the oscillations of the decay rates. For larger 
times, i.e., t 3> A^^, it increases in time as the sum of the two exponentials, with 
rates 7(^0 — fi, 00) and ^{uiq -\- Q., 00) respectively. It is easy to see that the smaller 
the value of A in the spectrum, the smaller the stationary decay rate 7(0^0 + 00), 
while "f{ijjQ — ri, co) — a does not change with A. A consequence of this point is the 
possibility of having situations where, starting with the atom initially excited and 
the cavity with zero photons, while the state has completely decayed, the 

population of keeps a constant value for long times, namely the ones in the 

interval ^{luq — Q,,oo)^^ <C t <C 7(^1^0 + fi, 00). A consequence of this fact is that the 
population of the atom is trapped in the excited state for an amount close to 25%. 

One may wonder how strongly this effect could depend on the particular model 
of spectrum we have chosen. In fact, the possibility of having a certain amount 
of population trapping in the system depends on the stationary decay rates only, 
so that the effect of population trapping involves only the values of the spectrum 
at the Bohr frequencies of the transitions of interest, and not the shape of the 
spectrum over all the real axis. To clarify this point, in the next section we show 
the behavior of the same populations when the spectrum of the environment has a 
shape different from the one considered in this section. 
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4. Evolution of the system in the case of a structured spectrum 

As a second model for a non-flat spectrum, we take a simple model o f stru ctured 
reservoir, consisting in a Lorentzian background with a Lorentzian gaji^^^^ 



1 



aiXl 



1 



(9) 



27r {uj\ — ujY + X\ 27r {uj2 — uj ) -r ^2 

where the gap is given by the inverted Lorentzian peak at frequency u> = 0)2- 
From the calculation of the previous section, it is straightforward to show that each 
Lorentzian peak gives a similar contribution, but with opposite sign, so that the 
time-dependent decay rate of a transition of Bohr frequency to is given by: 



{lui - LU)^ + Xl 



1 



Al 



■ sin(u;i — uj)t — cos(u;i — uj)t 



-Alt 



a2A 



{UJ2 - W)2 + A2 



1 



UJ2 



Xo 



■ sin(ct;2 — i^)t — cos{lu2 — uj)t 



-\2t 



(10) 



In the following we will take ui — 102 — t-^o + and ai > a2 and Ai > A2: 
in this way we assume that the Bohr frequency corresponding to the transition 
\Eq) corresponds to the minimum of the spectral density, in analogy with 
what done in the case of a single Lorentzian peak. The choice of ai — 0.1 * 2J7, 
a2 = 0.099 * 2n, Al = 100 * 2n and A2 = 0.1 * 2fl gives a ratio 1/100 between the 
stationary values of the two decay rates of interest 7(0^0 i i7, cx)), a situation which 
is close to the ideal case of perfect population trapping. 



(a) 



7 0.1 



0.06 
0.04 
0.02 



(b) 




Fig. 1. (a) Time-dependent decay rates of (solid line) and (dashed line) as a 

function of t (in units of (2f2)^^). (b) Population of the atomic ground state as a function of t (in 
units of (2n)-i). 



As we see from Fig. l-(a), the short-time behavior of the rates is rather different 
from the case of the single Lorentzian peak, indeed in the case of the gap the decay 
rate of |i?i,+) initially increases in the same way the rate of |-Ei,-) does, and, after 
reaching its maximum value close to 0.1*2f2, it decays exponentially to its stationary 
value 0.01 * 2il. Anyway this difference in the behavior of the rates does not lead 
to observable effects in the dynamics of the atomic populations. Indeed, comparing 
the population of the atomic ground state in Fig. l-(b) with the corresponding 
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population in Ref. F9l we see that the behavior predicted by both the single-peak 
Lorentzian model and the Lorentzian gap model is exactly the same: for a long time 
a population trapping occurs in the excited atomic state, for an amount of about 
25%, no matter which of the two models one chooses. 

Summarizing, in the situation we have analyzed, the essential point to properly 
choose the model is to check the appropriate ratio between the stationary values 
of the decay rates. In this sense, the most important property is the value of the 
spectrum at the Bohr frequencies of the transitions of interest. On the other hand, 
the particular shape of the spectrum one chooses affects only the details of the 
short-time dynamics, especially the short-time behavior of the decay rates, but the 
main aspects of the dynamics of the atom-cavity system are not affected by its 
choice. 

We feel that these conclusions are quite general and that they apply to a wide 
variety of situations involving lossy cavities with non-flat spectra, as in the most 
advanced quantum electrodynamic systemi^. 
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